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5^ ' Abstract. 

jy« . Total scattering neutron powder diffraction measurements were performed on the 

tetragonal phase (a=6. 4202(5) A, c=9. 5762(12) A) of CI4. The experiments were 

followed by Reverse Monte Carlo (for POWder diffraction (RMCPOW)) modeling. 

C^ . Detailed analyses of the resulting particle configurations revealed that the observed 

diffuse scattering originates from the libration of the molecules. By examining the 

I ■ partial radial distribution functions a distinct carbon-iodine peak at 4.5 A is found, 

r^ , which appears as a consequence of corner-to-face mutual alignment of two molecules. 

O ' The occurrence of edge-to-edge alignments is also significant within the first carbon- 

I ^1 . carbon coordination shell. 

> 

Cn ■ PACS numbers: 61.05.fm, 61.43.Bn 
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1. Introduction 

Among the crystalline CX4 halides that do not contain iodine (CF4, CCI4, CBr^) 



U . monoclinic (CF4 P El [H I3 [5] , CCk^ and CBnff\), face-centered cubic {CCU^ and 

- - CBr^^) and rhombohedral phases (CF4[1], CC/4[6j and Ci?r/f[TO]) are common. Over 

a long period of time, the structure of carbon tetraiodide has been debated; initially 
it was classified as cubic (see the reference of [H]) or monoclinic [12]. analogous with 
the above carbon-tetrahalides. The crystal structure could eventually be solved in 1981, 
when it was classified as tetragonal (J42m), on the basis of single crystal x-ray diffraction 
data pjj. The sample in that work was synthesized by the reaction of ethyl iodide and 
carbon tetrachloride, using the synthesis of McArthur et al. |13j . 

Considering existing knowledge on XI4 type crystals [HI UHl [161 [El [13 [IE] that 
possess a cubic symmetry, this unusual structure has been explained by the close packing 
of / atoms and by an intermolecular I — I bond[rT]. The observed large value of 
the ionization potential, found by ultraviolet spectroscopy, was explained also by the 
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presence of this bond[2n]- Recently, this system was exploited as a model for packing 
tetrahedra[T9]. 

Another interesting feature of this class of materials is the possible presence 
of plastic crystalline phases, appearing as the rhombohedral phase of carbon 
tetrafluoridejl], and the face centered cubic phases of carbon tetrachloride [HI [21] and 
carbon tetrabromide[9j. The common feature of these phases is that while the centres 
of the molecules maintain translational (i.e., crystalline) symmetry, the molecules can 
rotate (for CF4, see [22]). This rotation is not free, but hindered by neighbouring 
atoms. As a result, short-range order exists between ligands belonging to neighbouring 
molecules. This manifests as diffuse scattering that becomes structured in reciprocal 
space, due to its relation to the average structure (e.g., for CBr^^, see [231 l2ll 125]). 

Detailed analyses were presented some years ago for some of the crystalline CX4 
compounds, applying Reverse Monte Carlo modeling [261 EH] and molecular dynamics 
simulation [27]. These works showed that correlations between the 'ligands' (CI atoms 
in CC/4 and Br atoms in CBr^) have a close relationship with those found for the 
corresponding liquid phase. 

During our systematic studies of crystalline phases of carbon-tetrahalides, diffuse 
scattering in carbon tetraiodide has been found; this initiated the present study. 
Since the analysis of the Bragg-pattern has already been done by Pohl[TT], the 
present investigation focuses on pair correlations between atoms and on orientational 
correlations; concerning these issues, standard crystallographic analysis does not provide 
complete information on the structural details of the system. 

2. Experiment 

The sample has been provided by Aldrich (97% purity) and was further powdered before 
the experiment so that random orientation of the crystallites could be guaranteed. 

Room temperature carbon tetraiodide has been measured by neutron powder 
diffraction on the former SLAD diffractometerj28] (Studsvik NFL, Sweeden). Using 
the standard wavelength of the instrument (1.119 A), the diffraction pattern of the 
sample could be recorded in the 0.29-10.55 A~^ range. Scattered intensities from the 
8 mm diameter vanadium sample container, background and a solid vanadium rod have 
also been determined. Standard data normalization and correction procedures [29] were 
performed by the CORRECT[30] computer program; absorption, multiple and inelastic 
scattering contributions have been taken into account during the process. 

As it turned out later, the sealing of the sample holder has leaked during the 
measurement, and a small amount of sample has been evaporated from the can. In the 
pattern several additional Bragg-peaks appeared, which are supposedly originate from 
the orthorhombic phase of I2 [311 [32] (spontaneous brake-up of CI4 molecules, due to 
the tight packing of large iodine atoms, is common). The most intense such 'alien' Bragg 
peak can be observed at 1.7A^^. 

The lattice parameters have been re-determined by the 'Fox' software [33j: they 



The origin of diffuse scattering in crystalline carbon tetraiodide 3 

were found to be a=6. 4202(5) A, c=9. 5762(12) A. The amount of the I2 impurity has 
been estimated by a joint fit of both compounds using the estabhshed structures [3T| 132] : 
the amount was found to be about 2.8%. This corresponds well with the guaranteed 
purity of the commercially available compound. This indicates that the scattering 
from impurities originates from the sample, and not from the accidentally deposited, 
previously evaporated vapour from the container wall. 

3. Structural modeling 

Computer modeling has been performed by using the Reverse Monte Carlo for 
POWder diffraction (RMCPOW) software[3l]. RMCPOW is capable of simulating 
both Bragg- and diffuse scattering in Q-space. The molecule has the shape of an 
ideal tetrahedron, with 2.157 A bond length between carbon and iodine atoms [35l l36] . 
The molecular structure was maintained during the simulation by 'fixed neighbours 
constraints' (fnc)[37]. Carbon-iodine and iodine- iodine distances within the molecules 
were constrained to the distance regions between 2.1 to 2.2 A and between 3.42 to 3.62 A 
, respectively. These constraints, apart from keeping molecules together, allowed them 
to rotate by series of single atomic moves. The intermolecular closest approach distances 
were set at 5.7 A, 4.0 A and 3.3 A for carbon-carbon, carbon-iodine and iodine-iodine 
pairs, respectively. 

Initially, a supercell of 10'^ times of the Bravais-cell has been constructed, which 
corresponds to 2000 molecules. Atomic positions in the unitcell were set as published in 
[TT] . using the previously determined lattice parameters. At the end of the simulation, a 
goodness-of-fit R^p = 9.54% has been reached and 20 independent configurations have 
been collected for further analyses. (Between two 'independent particle configurations', 
each atom in the simulation box has moved successfully at least once.) 

4. Results and discussion 

Diffuse scattering from a powdered crystalline sample may result from effects other than 
positional disorder in the bulk crystalline single phase. Here it is assumed that positional 
disorder is the only reason for the observed diffuse scattering; before discussing results 
from our calculations it is therefore necessary to consider some of the other possibilities. 
First, one might argue that the origin is the molecular iodine (I2) that is present 
in the sample as impurity; this is not really probable, due to its small contribution and 
due to the lack of an intense diffuse scattering in the pure phase measurements [31], [32] . 
Another hypothesis is that diffuse scattering might be related to the (small) size of the 
grains. However, no broadening of the Bragg-peaks has been observed here that could 
be related to small grains. As yet another possible reason, the presence of an amorphous 
phase may be assumed. Amorphisation has been reported for CCI4 (see, for example, 
[38]). which took place during vapor deposition to a cooled substrate. For Snl^ the 
presence of amorphous phase is reported at pressures over 25 GPa [39]. In contrast to 
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these (somewhat extreme) cases, the present study was performed at ambient conditions 
(and on commercial samples) and no trace of any additional amorphous phase could be 
observed. We therefore can conclude that the observed diffuse scattering is related to 
bulk phase structural properties of crystalhne C/4. 

At the start of the simulation, we tried to remove the Bragg-scattering part of 
crystalline I2 traces from the full pattern, by extrapolating the components of the Debye- 
Waller factor from the low-temperature measured data|32]. This resulted in a slightly 
better agreement with experimental data; still, some of the discrepancies remained, 
most probably due to the not-well-known magnitude of Debye- Waller factors. This is 
why we decided not to bias results with this complication and started the simulations 
without subtracting this kind of background. At the final stage of the RMCPOW 
calculation, an agreement shown in figure [1] was reached, which (with the R-factor value 
less than 10 %) may justifiably be called satisfactory. A possible reason for the remaining 
differences between measured and calculated data is the presence of Bragg-peaks of the 
small amount of I2] they are most apparent between 1.5 and 2.5 A~^ where they caused 
small modulations of the diffuse scattering part of the simulated pattern. 

The coherent neutron scattering lengths of carbon and iodine are 6.646 and 5.28 
fm[l0], respectively. Thus the contributions of partial structure factors to the measured 
total scattering structure factor are 5.7%, 36.4%, 57.9% for carbon-carbon, carbon- 
iodine and iodine-iodine pairs, respectively. C-I and I-I correlations may therefore be 
determined well, whereas the position of C atoms (and hence, C-C correlations) is 
influenced mostly by the applied fixed neighbour constraints. Since none of the 2a and 
8i Wyckoff sites restrict the reflection conditions (in comparison with the general 16j 
WyckofF site), the above contributions of the partials are valid for all Bragg-refiections, 
too. 

Figure [2] shows a condensed view of the supercell as reduced to one Bravais-cell. 
This picture represents not just the average structure but also the effect of thermal 
displacements, which latter manifests as the difference between instantaneous and 
average positions. Intuition based on inspecting the spreads of instantaneous atomic 
positions around their equilibrium values suggests that the major contribution to the 
measured diffuse scattering must come from the sizeable displacements of I atoms or, in 
other words, the libration of the molecules. 

4.I. Partial radial distribution functions 

Partial radial distribution functions (prdf), as calculated from the final configuration, 
are shown in figure [3l Among the carbon tetrahalides, the closest analogue of the 
ordered phase of carbon tetraiodide is the ordered phase of carbon tetrabromide; for 
the sake of comparison, rescaled prdf-s[25] of CBr^ in phase II (monoclinic) are also 
shown. Since the two structures do not belong to the same spacegroup, exact rescaling 
(connected to the different sizes of Br and I atoms) is not possible. The scaling factor 
along the 'x' axis was chosen to be 1.06 for carbon-carbon and for carbon-iodine pairs. 
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taking into account approximate positions of maxima of the carbon-carbon prdf. For 
iodine-iodine correlations a factor of 1.1 was chosen, taking into account the minimum 
positions between the first and second maxima. 

Comparing the two phases, the carbon-carbon prdf of carbon tetraiodide looks much 
more ordered; the main origin of this difference is the different crystalline structure (and 
perhaps the fact that the Ci?r4 data were measured only about 20 K below the phase 
transition temperature is also an issue in this respect). The C-C maxima in CI4 can 
be found approximately at 6.5 A, 9.1 A, 11.2 A, 13.1 A (and larger distances, too); 
these values nearly correspond to the shortest lattice side (6.42 A), the half of the body 
diagonal (6.59A), the 'xy' diagonal (9.08 A), the second neighbour half body diagonal 
(11.22 A), and the full body diagonal (13.15 A) values, respectively. 

It may be appropriate to discuss here briefly the sharp (unphysical) 'cut-offs' of the 
C — C and / — / partials at the shortest distances. Usually, this feature may have two 
sources: (1) the higher than real intermolecular cut-off and (2) non-physical short period 
oscillations in Q-space. Taking into account that the iodine-iodine cut-off distance is 
lower than the intramolecular low limit, the observed phenomena must be result of that 
the program tried to fit the Bragg-peaks of J2. 

Comparing the centre-ligand and ligand-ligand prdf-s with the corresponding prdf- 
s of CBr^ phase II, C/4 again appears to be a more ordered system. Concerning the 
short-range behaviour, I-I maxima are at 4.05, 4.4 and 6.6 A, where the first and the last 
ones are distinct. Although the peak at 4.05A appears, that should possess a crucial role 
in explaining the formation of the tetragonal lattice [HI [20], it is not so well-defined in 
terms of pair-correlations. In contrast, the carbon-iodine peak at 4.5 A looks much more 
characteristic. For assigning the observed peaks to local orientations, see section 14.21 

4-2. Angular and orientational correlations 

Angular distributions (figure H]) C-C-C and C-I-C were investigated in the following way: 
atomic triples contributed to the statistics if the distance between carbon atoms were 
less than 7.5 A and less than 4.85 A between carbons and iodines. These values are 
upper limits of the first coordination shells of the corresponding prdf-s (see section H?T|) . 
For C-C-C triplets maxima at 60, 90, 120 and 180 degrees appear, which correspond 
well to the local symmetry of the crystal. However, the spread around the regular angles 
is not the same: the narrowest maximum is the 180 degree one, which shows chain-like 
ordering of the 'a' sides of the unitcell and through the body diagonal. Less ordered is 
the body diagonal-'a' side combination (at 60 and 120 degrees) and the broadest one is 
between the 'a' sides. 

For C-I-C correlations, a maximum is found around 165 degrees (-0.965 in cosines) 
and two less intense ones at 120 and 90 degrees. The first one is at a larger angle than the 
crystallographic result pT| (153.6 degrees (cosine: -0.896)), showing that iodine atoms 
prefer to be close to, but not exactly on, the carbon-carbon connecting line. Since the 
intensities of the remaining two maxima are much smaller (note the logarithmic scale in 
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the figure!), supposedly they are spurious intensities arising from C-I correlations that 
can be found at the high end of the distance range. This kind of correlation is then 
responsible for the observed peak of the C-I prdf at 4.5 A. It is instructive to note that 
if we calculate the average angle from the known maxima of the partial rdf-s (C-C: 6.5 
A, C-I: 4.5 A and the bondlength of 2.157 A), then we obtain almost the same result as 
the one derived from crystallography. That is, the average value is actually somewhat 
misleading here; clearly, the distribution is the proper characteristics of the structure. 

For describing orientational correlations between pairs of tetrahedral molecules as 
a geometrical unit, the classification scheme of ReyjH] has been used. This is based on 
the number of atoms of each of the two molecules found between two parallel planes, 
which planes contain the centres of the corresponding molecules and are perpendicular 
to the line connecting the two centres. Based on this consideration, six classes exist: 1:1 
(corner-corner), 1:2 (corner-edge), 2:2 (edge-edge), 1:3 (corner-face), 2:3 (edge-face), 3:3 
(face-face). The asymptotic limits of the probabilities of 1:1 and 3:3 are the same; this is 
also true for the 1:2 and 2:3 arrangements. These equivalences come from the fact that it 
should be 4 atoms on average between the two planes. In the case of tetrahedral liquids 
this resulted in an alternating behaviour between the 1:1 and 3:3, and between the 1:2 
and 2:3 orientation probabilities (c.f. Refs. [HI HJl HSj HH US]). Similar behaviour was 
observed in the plastic crystalline phases of the corresponding materials [25| Hi] . 

Figure [5] shows the probabilities of each class as a function of the distance between 
the centres of two molecules. A general observation is that only 2:2 and 1:3 arrangement 
play major roles in the system, others are marginal. It should be noted that there are 
distance ranges where the carbon-carbon prdf contains only very few atoms in a given 
distance bin (which is quite normal for a crystal)- i.e., the definition of any 'distribution' 
would be meaningless. This is the reason why between 7.5 and 8.2 A no intensities are 
shown in figure [5l 

Edge-edge arrangements are the most probable ones: we can find them in great 
numbers around the maxima of the carbon-carbon prdf. However, the first maximum of 
the 2:2 distribution appears at shorter distances than the position of the first maximum 
of the C-C prdf. Comparing with the monoclinic phase of CBr^^, the 2:2 arrangement 
there is much less preferred at shorter distances; instead, the phase change from the 
ordered to the plastic crystalline phase of CBr^ is explained by the transformation of 
several classes into 2:2 in the plastic crystalline phase[25]. The large number of 2:2 pairs 
at the 'contact' (very short) distance seems to be a distinct feature of tetragonal CI4. At 
medium range (above 8 A), the two crystalline materials follow similar characteristics. 

Similarly, the Apollo-like 1:3 arrangements appear in the first shell with quite a 
large probability; in contrast to the 2:2 class, 1:3 pairs are most frequent around the 
exact position of the first maximum of the C-C prdf. According to the angular analysis 
above, the iodine close to the line connecting neighbouring carbon atoms would suggest 
a longer centre-centre distance (6.6 A) than it is observed (6.5 A): this kind of distortion 
is shown in figure El As a consequence of the crystalline symmetry, the 1:3 arrangement 
becomes very probable also at 13. 3A, which corresponds to the distance of the body 
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diagonal of the Bravais-cell. (The shght difference between observed expected distances 
seems to propagate along the body diagonal.) We found similar characteristics in the 
monoclinic phase of CBr^ at medium range. A general observation concerning the 
two most important orientations: along the body diagonal, 1:3 arrangements dominate 
whereas along the 'a' vector, 2:2 pairs are more abundant. 

Although their role is marginal, it is still worth spending a few words on the 
remaining four classes of orientations. 1:1 and 3:3, and 1:2 and 2:3 arrangements are in 
phase (not like in the corresponding liquid phases [13] ). which can only be explained by a 
well defined crystalline structure. These pairs appear with relatively larger probabilities 
where the intensity of the carbon-carbon prdf is low. The only exception is the 2:3 
arrangement that shows up with a considerable probability in the region of the closest 
molecule-molecule distance: their presence may be taken as deviation from the (nearly 
uniform) 2:2 alignment and thus, they may contribute to the observed diffuse scattering 
intensity. 

5. Conclusions 

We performed total neutron scattering diffraction measurements on the tetragonal 
crystalline phase of carbon tetraiodide. The experimental total scattering structure 
factor was modeled successfully by the RMCPOW|34j algorithm. As seen 
unambiguously from displaying atomic positions (see Figure E]), the observed diffuse 
scattering comes from an enhanced positional disorder of iodine atoms, i.e. from 
molecular librations. The previously suggested |11] intermolecular iodine- iodine 
correlation at 4.05 A has been found, although its manifestation in terms of the prdf 
is not very distinct. The carbon-iodine correlation centered around 4.5 A has been 
found to be significant; this feature of the corresponding prdf results from 'Apollo- 
like' (1:3) mutual arrangements. The cosine distribution of the C-l-C intermolecular 
angles suggests that iodine atom positions close to, but not on, the intermolecular 
carbon-carbon connecting line are preferred. Analyses of mutual orientations of two 
molecules showed that edge-edge and corner-face (Apollo-like) near neighbours are the 
most important and also, that these types of arrangements remain distinct at larger 
intermolecular distances, too. This behavior is completely different from that found for 
the monoclinic phase of C.Br/i[25]. 
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Figure 1. Experimental and simulated powder diffraction patterns of tetragonal CI4 
at 298 K. Circles: measured differential cross-section; solid line: RMC calculated total 
scattering intensities; light red line: RMC calculated diffuse scattering intensities. 
The red arrow marks the most intense parasitic Bragg-peak, supposedly from iodine 
impurities. 
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Figure 2. Condensed view of tlie simulated configuration (fO'^ unit cells), as projected 
into one Bravais lattice. Balls represent individual atomic positions. Black: carbon; 
magenta: iodine. The considerable spread of the iodine positions can be assigned as 
the major reason behind the observed diffuse scattering intensity. 
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Figure 3. Partial radial distribution functions of crystalline CI4 calculated from the 
simulated configuration. Upper panel: carbon-carbon; middle panel: carbon-iodine; 
lower panel: iodine-iodine correlations. Red solid line: CI4 RMCPOW intermolecular 
part; red dashed line: CI4 RMCPOW intramolecular part; light blue solid line: the 
corresponding RMC simulated j25j CBr4 phase II partials (the 'x' axes for these curves 
have been re-scaled, via multiplying by 1.06 (upper and middle panels) and 1.1 (lower 
panel).) 
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Figure 4. Cosine distributions of some intermolecular angles. Upper panel: nearest 
neighbour C-C-C angle. Lower panel: C-I-C angle distribution. (Note the logarithmic 
'y' axes.) 



The origir 



14 




Figure 5. Rey's constructions [H] of orientational correlation probabilities as functions 
of the centre-centre distance. Upper panels [(a) and (b)], black solid lines: 2:2 [(a)] and 
1:3 [(b)] arrangements; red solid lines: the corresponding probabilities as calculated 
from the average structure; light violet: the corresponding probabilities for ordered 
crystalline CBr4^^ ('x' scale multiplied by 1.1). Lower panels [(c) and (d)], black solid 
lines: 1:2 [(c)] and 1:1 [(d)] arrangements; red solid lines with circles: the corresponding 
probabilities as calculated from the average structure; grey tone lines: 2:3 [(c)] and 
3:3 [(d)] arrangements; blue solid lines with crosses: the corresponding probabilities as 
calculated from the average structure; olive line: the centre-centre prdf (shifted along 
the 'y' s-xis and re-scaled for a better visibility). 



